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Abstract 

The possibility of treating boundary conditions in terms of a bilocal dynamical field is formal¬ 
ized in terms of a boundary action. This allows for a simple path-integral perturbation theory 
approach to physical effects such as radiation from a time-dependent boundary. The nature 
of the action which governs the dynamics of the bilocal field is investigated for a limited case 
(which includes the Robin boundary conditions). 


1 Introduction 


The role of boundary conditions in field theory has recently been of considerable research 
interest. This has been explored in great detail in the case of the Casimir effect, for a variety 
of geometries pQ. The question obviously goes beyond the Casimir effect, which only pertains 
to the vacuum-to-vacuum amplitude. The impact of different boundary conditions on higher 
point functions is clearly of importance in diffraction theory and other physical phenomena. 
The general theory of self-adjoint extensions mm also allows for boundary conditions with 
negative eigenvalues for the operators of interest such as the Laplace operator; these are related 
to edge excitations and can lead to interesting phenomena mm- 

From the mathematical point of view, boundary conditions are necessary to render the 
problem well-defined with the needed self-adjointness properties. But from a physical point of 
view, boundary conditions are idealizations of dynamics on the boundary. Thus the dynamics of 
material particles which constitute a metal plate would lead to nonzero electrical conductivity 
and, in an idealized limit, would give the standard Dirichlet and Neumann conditions used 
in Casimir calculations with the electromagnetic field. Thus it is physically meaningful to 
have boundary conditions which can vary from point to point on the boundary or which can 
evolve with time, in accordance with some dynamical principle. Effectively, one must treat 
boundary conditions as additional dynamical fields with their own action and time-evolution. 
In particular, one may ask how the dynamics of boundary conditions for a given set of fields is 
modified by the back reaction of the fields themselves. In this paper, we start exploring some 
of these questions. 

In considering diffractive contributions to the Casimir effect, we have recently developed a 
boundary action approach m-m- The strategy is to obtain a lower dimensional held theory 
dehned on the boundaries, as a functional of the boundary value of the held, by integrating 
over the bulk helds. This boundary action makes it straightforward to incorporate the effects 
due to edges and apertures on the boundary as part of the integration over the boundary 
values of the helds. The formalism was also extended to include general boundary conditions 
as allowed by the von Neumann theory of self-adjoint extensions [9]. In this case, an integral 
kernel in the boundary theory encodes the allowed boundary conditions. This formulation 
has the advantage of recasting the entire discussion of boundary conditions and how they can 
change with space or time as a lower dimensional held theory which lives on the boundary. In 
particular, one does not have to modify the mode expansions of the bulk helds as one changes 
the boundary conditions. Thus boundary conditions can indeed be treated as a dynamical 
held. In this paper, we will take the next logical step and consider how the bulk helds can 
affect the dynamics of the boundary conditions. 

The paper is organized as follows. In the next section, we will briehy review the bound¬ 
ary action. In section 3, we consider radiation from the boundary due to time-dependent 
Robin conditions. This problem is closely related to the dynamical Casimir effect since time- 
dependent boundary conditions can simulate moving mirror effects HO]-US]. The calculation 
of the radiation from boundaries with time-dependent Robin boundary conditions was done 
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by using Bogoliubov transformations to define a new set of creation and annihilation operators 
for the field mm- Our method will be much simpler, involving the perturbative expansion 
of a path-integral for boundary fields, which can easily accommodate arbitrary dimensions, 
higher order effects, corrections due to bulk interactions of fields, etc. 

In section 4, we consider the effective action for boundary conditions, i.e., for the integral 
kernel or the bilocal field on the boundary which encodes the boundary conditions, which 
obtained by integrating out the fields. This is a rather involved problem, our calculation is for 
a limited choice of the bilocal held. Nevertheless, it illustrates the method and some general 
features of the action for this held. A more complete analysis is clearly called for, this is 
currently under investigation. 

The paper concludes with a short discussion. 


2 General boundary conditions 


We review the boundary action briehy considering a free scalar held in a region V with 
the usual kinetic term of the form Jy(d(j))‘^/2. On the boundary, the combination i-p -|- idn^, 
where is the boundary value of cj) and dnP is its normal derivative, is to be viewed as an 
element of a Hilbert space of L^-functions. The general boundary condition, according to the 
von Neumann theory of self-adjoint extensions, is mm 


if + idnP = U{ip-idnp) 

{(f + idn<p){x) = (j) U{x,y) {(f - idn(p){y) 

Jv 


( 1 ) 


where 17 is a unitary operator on the boundary Hilbert space; this is made more explicit in 
the second line by writing U{x,y). We can also rewrite the boundary condition as 


dnP = -i 


U-1 
U + 1 


if — —1C <p 


(2) 


where /C is a hermitian operator; it corresponds to the Cayley transform of U. The limit 
/C —>■ 0, equivalent to t/ = 1, gives the Neumann boundary condition, while /C —)• oo gives the 
Dirichlet condition, as seen by dividing ([2]) by /C and taking the limit /C —)■ oo. The case of 
fC being a constant (proportional to the identity on the Hilbert space) is the Robin condition. 
These are special points on the space of boundary conditions; clearly, in general, there is a 
much larger class of choices. 

In [6], we considered a two-step evaluation of the Euclidean partition function. While 
variants of this, including the possibility of interacting fields were also considered in [7], for the 
present discussion, we will consider a free massless scalar field theory with a standard action 
^{4’) ~ ^ f (^4^)^- Further we consider a boundary which is a plane normal to the xi direction. 
We can then write the partition function as 


Z = 



(3) 
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with the field 4>{x) parametrized as 


4>{x) = / i^{x')n ■ dx'G{x',x) + r]{x) (4) 

J x' 

where G{x',x) is the Green’s function for the Laplacian in the bulk with Dirichlet boundary 
conditions, n is the unit vector normal to the boundary and (p is the boundary value of cf). The 
field r]{x) also obeys the Dirichlet condition ry = 0 on the boundary. The parametrization ([4|) 
ensures the appropriate boundary behavior for (j). Integrating out rj one obtains 

Z = det(-n)y j [dip] (5) 

where denotes the boundary action. It is of the form 

= ^ [ ^{x')M{x',y')ip{y') 

^ Jx',y^ 

M{x\y') = n- da;'n- dy>G{x',y') (6) 


The general boundary condition ([2]) can be incorporated by using the augmented boundary 
action [9] 

^ / ip{x') [M{x',y')+IC{x',y')] ip{y') (7) 

The parametrization (jlj) does not imply any boundary condition on (j) since the boundary value 
ip is unrestricted. One may even think of the plane as a fictitious plane where we have (j) = ip. 
The boundary conditions on cf) are imposed by the boundary action ([7|), when we integrate 
over ip, after the choice of /C. As said before, /C = 0 corresponds to Neumann and /C —)• oo to 
Dirichlet conditions. 

To see how action ([7]) comes about, we start by considering the normal derivative of the 
field. Consider formulating the functional integral for Z by discretizing the coordinate along 
the normal direction, say, the xi-direction. The Euclidean action which enters (l3|) then has 
the form 


5 = Sim) = ^ 


'3^T 


d^x 


{(pN - i4>N-l - (pN-2Y _^ 

XN — XN-l XN-i — XN-2 


( 8 ) 


where (pN = ip is the boundary value of (p and the superscript T indicates components tangential 
to the boundary. Integrating out the fields cpi^ keeping (pj\f = ip fixed, we get the boundary 
action Ssi'p) given as 


q-SbM 


N-1 

exp(-S'(y5,{0i}) 

1 


(9) 


If we functionally differentiate e ^ with respect to (/? = (pN, we find 


<5 _s ^ -S - (pN-l) -s 

- 6 = - C =-6 

5ip d(pN Xn — XN—1 


-dnip e ^ 


( 10 ) 
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In other words, the normal derivative can be obtained as the result of functional differentiation 
of the e~^ with respect to the boundary value of the field. This is a key result for us as it gives 
a way to express the normal derivative entirely in terms of the boundary action. 

We can now see that the general boundary condition ([2]) can be obtained for the remaining 
integration over ip if we use the augmented action Sb{'P,IC) from ([7]). For this consider the 
identity 


j [dip] exp{-SBi<p,IC)) {dn<p) = j [dip] exp j ipJCip- {dnpi) 


= I [dip] exp ipJCip 


] e-5B(¥>) 

6ip^ 


j [dip] exp 


= I [dip] exp ( -^ / ipJCip] {-K.ip)e 


( 11 ) 


where, in the last line, we have done a partial functional integration. The first and last steps 
of this equation show that 


J [dip] idnP> + JCip)=0 (12) 

This justifies our argument that the augmented boundary action ([7]) incorporates the general 
boundary condition ([2]) . 

It is worth pointing out that the use of the action ([7|) to take care of boundary conditions 
has some advantageous features. The field r] and the Green’s function G{x',x) in Q obey a 
fixed boundary condition, namely, the Dirichlet condition. The real genesis of various boundary 
conditions for the held ([) is transferred to ip and the action ([7]) which controls it. This makes it 
much easier to analyze change of boundary conditions, including time-dependence, dynamical 
determination of boundary conditions, how back-reaction from other helds can modify the 
boundary conditions, etc. We can treat the boundary conditions effectively as a boundary held. 
The fact that this approach is viable, namely that the boundary action ([7|) does reproduce the 
physical effects of different boundary conditions correctly, is clear from calculations of the 
Casimir energy pin]. We also point out that a similar idea of dynamically implementing the 
standard Robin boundary condition in terms of a boundary action was also suggested in [H], 
although our formalism is somewhat more general. 


Going back to (I7|) , we see that we can impose different conditions on different subsets of the 
boundary. Regions with Neumann condition will have JC = 0, while regions on the boundary 
where /C —)• oo must have ip = 0. We may regard the hermitian operator JC as an integral kernel, 
as already indicated in ([7]). Thus, generally, boundary conditions correspond to a bilocal field 
JC{x',y') on the boundary and the general form of the boundary action is 


Sb = SB{ip,JC) + Sfc (13) 

where Fyc is the action for the field lC{x',y'). We have obtained the form of SBi<p,JC), but 
the question arises: What is the form of the action Sjc for the bilocal field IC{x,y)7 A related 
question is: To what extent is it advantageous to think of IC{x,y) as a dynamical field? 
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We will take up the second question first. In the next section, we consider some examples 
which highlight the value of viewing /C as a held. We calculate the radiation from a time- 
dependent boundary condition. Such boundary conditions can be engineered in many ways: 
For example, for the electromagnetic held, one can consider a metal plate as the boundary and 
modulate the conductivity. This can also be viewed as the problem of moving mirrors where 
the effect of the movement can be transferred to time-dependence of the boundary condition 

m- 

Regarding the hrst question, clearly, the form of S)c is to be determined ultimately by the 
physics of the material of the boundary. However, we can gain some insight into the nature of 
Sic by calculating corrections to it due to the helds ip, i.e., by integrating out ip. This is familiar 
from standard held theory. If we integrate out one set of helds and obtain the corrections to 
the action for another set of helds, we can postulate that the action for the latter must at 
least have terms of the appropriate structure (as monomials of the helds and their derivatives) 
to provide proper renormalization for the terms generated by integrating out the hrst set of 
helds. A simple and age-old example is when we consider fermions coupled to a scalar held. 
Integrating out the fermions shows that we must have the usual kinetic term, mass term and an 
additional quartic self-coupling for the scalar for proper renormalization. This is the strategy 
we will follow here. We will consider the nature of the action for JC generated by integrating 
out the ip helds. 

Finally, although we will use a simple free held theory example for the calculations to 
follow, we may note that the formulation is quite general and can easily handle interacting 
theories. The construction of the boundary action for the general case was briehy discussed in 
[7]. The part of the boundary action without the KLip^ term is dehned by 

where r] vanishes on the boundary and (j) takes the value ip on the boundary, (j) is extended 
into the bulk from its value ip on the boundary in a unique way, so that there are no additional 
functional degrees of freedom associated with it. The result (flTl) can be phrased as follows. 
Consider the generator r[y] of the IPI vertices of the theory, which is defined by 

exp (-r[x]) = j [dr]] exp + v] + j (15) 

Here we will consider evaluating F such that the field r] vanishes on the boundary. This 
means that the Green’s functions occurring in, say, a perturbative expansion of F will obey the 
Dirichlet conditions. So far the field x is an arbitrary field. We now choose y as a particular 
solution of the equation of motion 

<5r , , 

— = 0 16 

5x 

obeying the condition y —>• (/j on the boundary. Comparing (fT4)l and (fT5]l with this condition, 
we see that the integral in (|15p gives the boundary action. In other words, 

Sb['p] = r[y], subject to — = 0, and y —)• on the boundary (17) 
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The full boundary action can then be written as 


1 


Sb = r[x] 

This gives the boundary action including the effect of interactions. 


+ - / (p{x)IC{x ,y)(p{y) + Sic 

{5r/5x)=o,x^v 


(18) 


3 Radiation from boundary 

We will consider a single flat boundary, say at xi = 0, with Robin boundary conditions for 
the field cj), with the parameter k taken to be a local boundary field. In other words, we take 
)C{x,y) = K,{x) d^^\x — y). We use a notation where the coordinates of a point in the bulk 
are denoted by (a;,a;i), where x corresponds to directions tangential to the boundary surface, 
including the Euclidean time direction and xi > 0 is along the direction perpendicular to 
the boundary. One immediate consequence of time-dependent boundary conditions is radia¬ 
tion from the boundary; this has been calculated for scalar fields and time-dependent Robin 
boundary conditions in (H-l) dimensions |11] by carrying out a Bogoliubov transformation on 
the fields. From our point of view, the Robin boundary condition corresponds to a boundary 
action term j KyP'l2 for a scalar field. Thus we may view the radiation as a decay process 
for the K field, k ^ (ftp. It is then straightforward to calculate this in simple perturbation 
theory. The physical observation is that we have two detectors placed somewhere in the bulk, 
say at points x and y, which show the transition corresponding to the absorption of two of 
the (/3-particles. The boundary fields propagate into the bulk as given in Q. Let J designate 
a source operator associated with the detector whose matrix element effects the transition 
detecting the (/3-particle. The amplitude for the process is then given by 

A = J Jfi{x,xi) JpiPy,yi) {(f){x,xi)(f){y,yi)) 

= J xi)Jf,i,{y, yi){(p{x')(p{y'))dvG{x', x'l, x, xi)U;=o dvG{y', y[;y, yi)|y;=o(I9) 

where {ip{x')(p{y')) is the propagator for ip as given by the boundary action and G{x, xi; y, xi) is 
the Dirichlet Green’s function. In order to derive the transition amplitude ()I9p corresponding 
to the absorption of two particles of energy uj and uj', it is necessary to analytically continue 
to Minkowski space. We briefly outline how this is done. For a single boundary, taken to be 
at xi = 0, the Euclidean Green’s function is given by 

df^k dki 


GE{x,xi-,y,yi) = 2 J J —j-^^^sm{kiXi)sm{kiyi) (20) 

The analytical continuation to Minkowski signature is done by using the substitutions 


XE = {Xd+l,X\\) ^ XM = {ixo,x\\) 

kE = {kd+i,k\\) ^ kM = {-ikoAw) 


( 21 ) 
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and the ie prescription. Here Xd+i is the Euclidean time and x\\ denotes the spatial directions 
tangential to the boundary. Using (|2T]) we find 


dx^G{x,xi\y,yi)\x^=o =-i / dy{k) 


dki 


1-00 kf + k'^^ — k^ — ie 
where d/r(/c) = dkod^\/{27r)'^. Using contour integration we find 
dx,G{x,xi-,y,yi)\x,=o = -i / dy,{k)e'‘ 


Jk(x-y) 


where 


Further 


ki = < 


if ^0 > ^\ 


i\lkn — Ajq if k?i > k^ 


( 22 ) 


(23) 


(24) 


M{x,y) = dx,dy,G{x,xr,y,yi)\x,,y,=o =-i J dy{k)e"’^^^ ^hki (25) 

where ki is given in (I24p . The Minkowski boundary action for the field ip is thus given by 

Sb{(P, ^) = -^ / ‘fix) [iki + ip{y) = -^ / [ih + + SK{x)]^y ip{y) (26) 


where we have separated out a constant term kq which is space-time independent. The last 
term in (I26p can be treated as the interaction part in a perturbation scheme. The propagator 
for (f is then given by 


{ip{x)ip{y)) = {ip{x)ip{y))+ I {ip{x)ip{z))^^{-i6K{z)){ip{z)ip{y)) 


K,q ^ 


' K,q 


+ 11 {ipix)(p{z))^^{-i6K{z)){ip{z)ip{z'))^^i-i6Kiz')){ipiz')(p{y))^^ + --i27) 

! z J z' 


where 


{ip{x)ip{y))= -i j Mk) 


^ik{x-y) 

Ko + iki 


(28) 


Using (EHI) and (I27P we can in principle calculate the transition amplitude in (fTUI) in ar¬ 
bitrary dimensions and to all orders in 5k. The calculation is particularly straightforward in 
(U-l) dimensions. In this case we find that the transition amplitude for the absorption of 
particles of energy lo, w' by the detectors is given by 


A = 


\^\)df'i'i+ 1 I) 

(ko -I- i\Lo\){Ko + i\ui'\) 


—5k{—oj — oj') + 


5k{—u} — — to') 


2tt 


(kq + 


-(-!)’ 


/ 


d^+'> 5 k{-UJ - e(^))5K(^W - ^(2)) . . . ^^(^(n) 


27r 


27r 


(Ko + ilC^^)|)---(/«o + *|C^”)|) 


(29) 


where J and 5k in (1291) refer to the appropriate Fourier transforms. In order to calculate the 
frequency spectrum of the particle number observed in the detectors we have to square the 













amplitude and sum over the final states. The last step is a bit tricky since we are summing 
over the final states of the detector. For a perfect detector, this process must be equivalent 
to integrating over the final particle phase space. This means that for the summation over a 
small range of final states we can use 

(30) 


A more detailed argument justifying (I30p is given at the end of this section. Using this and 
integrating over oj' we obtain the particle number distribution as observed in one detector as 


dN{ui) = duj 


2u f doj' 


a;' 0 (a;') 


TT 


2tt 

Sk{—u} — — u}') 


—6k{—u) — to') + j 

.-i-irj 


2-k 


{kq + 


d^W 6k{-uj - - ^( 2 )) • • • - w ') 


27r 


= duj 


2 a; r duo' uo'Q{uo') , , , 

TT J 2-k {kI + Uj‘^){kI + ^ ^ \ 


2 


(31) 

(32) 


Expression (13 ip is in agreement with previous results for particle production rate from time- 
dependent Robin boundary conditions. In m the frequency resolved particle production rate 
was calculated using Bogoliubov transformation for the fields up to second order in 6k. This 
agrees with ((3Tp with the notational identification, k = — 1 / 7 , kq = — 1 / 70 , 6k = —6-^/^^. 
Some higher order corrections, argued to be relevant in experimental setups for detecting 
dynamical Casimir effect m, were calculated in m, again using Bogoliubov transformations. 
Our expression ()3ip is in agreement with these results as well, after the appropriate translation 
in the notation. Finally, we note that the term in the square brackets in (131 p is the T-matrix 
if one regards 6k as a potential energy term for the ip on the boundary. This is made explicit 
in ((32p . 

The approach we used in calculating the radiation intensity is quite general and is applicable 
in arbitrary dimensions. We will now outline the corresponding derivation in (3+1) dimensions 
to order (d/t)^. In this case 6K{x,t) is in general a function of both space and time. Following 
the same procedure as before, we find that the transition amplitude for an absorption of two 
particles of frequency uo and uo' by the detectors when they are placed far away from the 
boundary is 


A = -i Jfi{LO,k\\,ki)Jf'i’{‘^',k'\\,ki) 


6k{—uo — uo', —{k + A:')||) 
(ko + ijkiDlKo + i\k']\) 


+ ••• 


(33) 


where ki = ya;^ — and /c( = y — k'\y This is the far-held contribution. If the detec¬ 
tors are placed very near the boundary, there will be a near-held contribution. The latter 
corresponds to the contribution from the range of k\\ such that > uo'^. We are ignoring this 
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for detectors placed far from the boundary. (The distinction between what is far and what is 
near is controlled by the wave length of the particle.) To continue along the lines of (1331) . the 
summation over the final states gives 


J{u,k)J*{uj,k) 

f 


2 Wfc 


d^k 

(27r)3 


(34) 


This is the higher dimensional analogue of ([30]l . Squaring the amplitude ([331) . using 


and 


integrating over u' with k[ = y — k'u we find 


dN (io) = 2uj 


d?k 

(27r)3 



J'^duj' \6 k{-uj - u',-{k + 9)^)^^ 
co"^ — P || (^0 ~ (^0 + — k'^^) 


(^ 


(35) 


If 6K{t) depends only on time and is space independent, (|35]) can be simplified further. We 
find 


dN{u}) = A 2uj 


d?k io'^duj' 


\5k{—u} — uj')\ 


l \\2 


( 27 r )^ _ ^2 (^2 _|_ ^2 _ fc 2 )( ft ;^ 2 _ ^ 2 ^ 

where A is the area of the boundary. 


(36) 


The basic idea in demonstrating the calculations (l3T]) and (l35l [SG)) was to phrase everything 
as a standard perturbative calculation in field theory, so that it brings out the nature of k as a 
boundary field; it also shows how one can easily include corrections, due to back-reaction from 
the bulk fields as well as multiple emission processes. 


Before we close this section, for completeness, we give the argument justifying (l30]) and (IMl) . 
These are really a standard result although not expressed in this form in most calculations. 
The amplitude for the emission (or absorption or scattering) of particles is usually written in 
the form of the integral over all coordinates of the product of a vertex function V{xi,X 2 • • •) 
and the single particle wave functions. 


A 



V{xi,X 2 , • • • 


,Xn,x)Uk, {Xl) Uk^{x2) ■ ■ ■ 


^ik-x 

y/2LOkV 



^ik-x 

y/2uJkV 


(37) 


where rt^’s are the wave functions (or their conjugates as needed), and in the second expression 
we have abbreviated the integral since we want to focus on one emitted particle of momentum 
k. The field is taken to be free and enclosed in a spatial cubical box of volume V with, say, 
periodic boundary conditions (with 17 —)• oo eventually). The summation over final states in 
the square of this amplitude gives the factor V{d^k/{2 tt)^), so that the phase space measure 
from the particle under consideration is 


1 d^k 1 d^k 

2ojkV (27r)3 2ojk (27r)3 


(38) 


There is another way to think about this process. The particles are emitted from the 
interaction region and absorbed by some detectors. One can consider this whole process; in 


10 



















fact, this is even more physicai as this is exactiy what is done in any experimentai observation. 
The emitted particie propagates from the boundary to the detector (taken to be at spacetime 
point y, with > x^) and so the ampiitude for this is 


A 


F{x)G{x,y) J{y) 


'x,y 



J{uj, k) 



^ipa{xo-yo)-ip-{x-y)) 

‘2ujp 


J{(X, fc) ^iijyO-ik-y 
27r 


(39) 


Upon squaring the ampiitude (1391) and summing over a smaii range of hnai states (which are 
now the hnai states of the detector), we get a factor | Jp/4a;^. This factor wiii be sensitive 

to the detection efficiency via the | Jp, but we can consider a perfect detector where this shouid 
give the same resuit as we obtained for the free particie. The agreement of this factor with 
(|38l) then shows that we need (IM|) for a perfect detector. 

This second way of thinking about the process has some advantage for us. For the hrst 
method, we wiii need wave functions with correct boundary conditions with JC inciuded. Whiie 
this is not particuiariy difficuit, our approach has been to use the boundary action directiy, 
so it is interesting to bypass the wave functions. The formuia (I34p re-expresses the needed 
summation as a property of the detector, and so the same resuit can be used for our case as weii. 
The same detector can be used for the particies emitted by the boundary. (Strictiy speaking 
there couid be a smaii back-reaction effect on the detector when a boundary is introduced, 
but to the order we are interested in this is not important.) As for the propagation of the 
particie from the boundary to the detector, we aiready have that taken care of in terms of 
dvG{x',x'^;y,yi). 


4 The effective action 

We now turn to the action Sjc for the biiocai boundary fieid JC{x, y). As expiained eariier, our 
strategy is to integrate out the (/9-heid and identify the kind of terms which wouid be generated. 
We then expect that such terms shouid exist in the action for JC. The fuii investigation taking 
account of the biiocai nature of JC wouid be rather invoived. Here we make a first attempt by 
considering /C’s of the form JC{x, y) = k,{x) S{x — y) = kq 6{x — y) + 6k{x) 6{x — y). With this 
choice of IC{x,y), we can easiiy integrate out the fieids cp in the boundary action and obtain 
the effective action for k as 

A5eff = ^ Tr iog(M + acq + Sk) (40) 

with M = '/IA. Here we have reverted to the Euciidean signature; since we are embarking 
on ioop caicuiations, it is easier to do this. We are aiso considering the case of the entire 
d-dimensionai boundary contributing; we do not assume Dirichiet condition on any part of the 
boundary. (Such regions, because we need k oo, cannot be treated in a perturbative fashion. 
We will need to use a modified M with expanded in mode functions with support only in 


11 





non-Dirichlet regions.) Also, to minimize clutter in equations, we will use the vector notation 
for momenta and coordinates; they will designate components tangential to the boundary. 

Expanding (jdni) in powers of 6k, we get the first correction as 

~ 2] “ ^(2wy 

The part of 6k which is independent of x can be absorbed into kq, so that, without loss of 
generality, we can take f d‘^x6K{x) = 0. Thus, we can take = 0. 

The term which is quadratic in 6k is given by 


Jk(x-y) 




+ Kq 


6k{x) 


(41) 


- -- 


/ 


d<^p 


d^xd^^y 


6k{x) 


^ik{x-y) ^ip(y-x) 


II 


d!^xd°'y 6 k{x 


(2n)H2^r ' 

dH 


6K{y) 


.( 2 vr) 


^-il{x-y) 


fd{l) 


6K{y) 


(42) 


where 


fd{f) 


I r d^p 1 

{\J{p + iY +^o) {^/W + ^^o) 


(43) 


The function fdi}) is even in I and we will see it has an expansion in terms of P, where I = \l\. 
Below, we evaluate ([451) in the case d = 3, 2,1. Carrying out the angular integrations in (H5|l . 
we get 


fd=3{l) 


1 f°° dpp 

8vr2/ Jq {p + Ko) 


l + p-\l-p\ 


Kq log 


/ Ko + l+p \ 
\Ko + \l-p\) 


(44) 


The above integral is divergent. This is due to the fact that fd{l = 0) in (f43l) diverges for 
d >2. After introducing an upper cutoff A for the p integration we find 


fd=3{l) 


1 

Svr^/ 



dpp 

(p + Ko) 


2p 


kq log 


/ Kq + / + p \ 
\kq + 1-p) 


+ 


j: 


dpp 


{p + Ko) 


21 — Kq log 


Ko + l + p 

Ko+p-l 


(45) 


For the effective action, we are interested in a small 1/Ko-expansion. 
evaluated in a straightforward manner to obtain 


h-M = 


2 A + 2ko — 4ko log 


A + Kq 
K o 


+ 


1 


72tt‘^ko 


The integrals can be 

■P +••• (46) 


where the ellipsis indicates terms of higher order in powers of P. Using this expression for 
in (I42p . we see that the first term of (1461) is like a mass renormalization, while the second 
term is like a wave function renormalization. Since the mass renormalization is divergent as 
A —oo, we have to postulate that the starting action has a term \ j p (h/i)^. Even though it 
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is not forced by the divergence structure, we can also consider adding a term proportional to 
f (VSk)^. Thus we consider the starting action Sofc 

Soic = Yj ( 47 ) 

With the addition of the contribution from (|42p . we get the effective action 


Z 

'2 


^ / 5k( — V‘‘ + li,ren)SK 


l^ren — 


Stt^ 


Z — Zq + 


1 


727 r 2 Ko 


2A + 2/to ~ 4/to log 

+ ••• 


A + /to 

/to 


+ 


(48) 


In (I48p . we could consider eliminating Z by scaling Sk to get a canonically normalized kinetic 
term for 5k, but this will redefine terms with higher powers of 5k in the expansion of MOj) as 
well as the | / 5k term. So we will leave the action ([18|) as it is. 

In the case of d = 2 we find 




log 


A + /to 

/to 


- I 


+ 


DGvr/tQ 


+••• 


This implies that the corresponding renormalized parameters for d = 2 would be 


Pren — A 


1 

dvr 


log 


A + /to 

/to 


- 1 


+ 


(49) 


(50) 


Z — Zq + 


Similarly in the case of d = 1 we find 


DGvr/tQ 


+ 


There are no divergences in this case and the corresponding renormalized parameters are 


(51) 


Mren — A 


T ■ ■ ■ 5 Z — Zo + — — 5 - + 


27r/to 


127r/tg 


(52) 


The calculations we have done are just the beginning in elucidating the nature of the action 
for JC. Even with the form of JC we have chosen, there are higher order corrections possible. 
Further the full bilocal nature of JC will bring in further complications and new features. Clearly 
continued investigations into these questions are needed, which we propose to take up in future. 
For now, some observations on the nature of K might be useful. 

The topology of the unitary transformation in ([T]) is an interesting question. It is a well 
known result, namely Kuiper’s theorem, that all the homotopy groups of the set of unitary 
transformations on an infinite dimensional Hilbert space % vanish. Therefore, a priori, it would 
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seem that we have only trivial topology. However, there are operators U which differ from the 
identity only on a finite-dimensional subspace of the Hilbert space 71, or phrased differently, 
of the form U = 1-1-0, where 0 is a compact operator. The homotopy groups of such unitary 
operators can be nontrivial in general. The question of which type of U's will be relevant for 
a given physical situation has to be answered on physical grounds such as requiring finiteness 
of energy. In other words, it will be determined by Sx.. If we consider restricted classes of t/’s, 
we can certainly get nontrivial topology; this has the potential to lead to interesting solitonic 
excitations of the boundary field. For such cases, staying within the chosen class, one may not 
be able to use simple boundary conditions such as the Dirichlet or Neumann ones as topology 
forces variation of JC{x,y) with x,y in a certain way. Clearly, the calculation of the action is 
crucial to discuss such questions. The case we have considered for the calculation of the action, 
namely, IC{x, y) = k(x) 5{x — y) does not quite fit into the case of C/’s of the form U = 1 -|- 0, 
with 0 compact. So new techniques are needed for evaluating the action in the general case. 

5 Discussion 

Boundary conditions can in general be encoded as an integral kernel. This kernel can be viewed 
as a bilocal dynamical field on the boundary. Physical effects due to boundary conditions, such 
as radiation from the boundary for time-dependent boundary conditions, can be calculated 
using standard field theory techniques such as perturbation theory. This is illustrated in 
section 3. The nature of the action which governs this bilocal field was considered in section 
4. By integrating out the bulk fields, we can obtain the general form of the action. This is 
worked out for a limited class of boundary conditions (which included the Robin case). The 
more general question of the action governing the bilocal field is under investigation. 

There are also more general situations than what we have considered here to which this 
method can evidently be applied. Moving boundaries, multiple boundaries, temperature de¬ 
pendence are some obvious examples. Situations which require nonlocal boundary conditions 
for which the bifocal nature of K,{x, y) and the action governing it are fully operational would 
be another important case. These will be taken up in future work. 

We thank Dan Rabat for a critical reading of the manuscript. This research was supported 
in part by the U.S. National Science Foundation grants PHY-1213380, PHY-1417562 and by 
PSC-CUNY awards. 
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